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Abstract 

We present a nonperturbative QCD calculation of diffractive vector meson 
production in virtual photon nucleon scattering at high energy. We use the 
nonperturbative model of the stochastic QCD vacuum which yields linear 
confinement and makes specific predictions for the dependence of high-energy 
scattering cross sections on the hadron size. Using light cone wave functions of 
the photon and vector mesons, we calculate electroproduction cross sections 
for /3, (j) and J /ip. We emphasize the behavior of specific observables such 
as the ratio of longitudinal to transverse production cross section and the 
t-dependence of the differential cross section. 
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I. INTRODUCTION 



Exclusive vector meson production by real and virtual photons is a good probe to inves- 
tigate the physics of diffractive scattering. Whenever the coherence length of the photon is 
larger than the proton radius, it is preferable to study the process in the proton rest frame or 
in the center of mass frame where the virtual photon can be considered as a hadronic system 
composed of partons. In this case the photon-hadron interaction is then in many respects 
similar to hadron-hadron collision. In addition, it offers the possibility to vary the polariza- 
tion and virtuality of the photon and thereby manipulate the light cone wave function of 
the incoming state. "The experimentalist can make hadrons of arbitrary sizes." 

In our approach, we attack the problem as a genuine nonpertubative one. We use the 
model of the stochastic vacuum which has been adapted to high-energy hadron-hadron 
scattering in Ref. 0], applying the general scheme developped by Nachtmann for the 
separation of the large energy scale from the small scale of momentum transfer. The model 
of the stochastic vacuum gives satisfactory results both in low and high energy physics. It 
yields a rather simple geometrical picture for a single gluonic flux tube 0]. The same mech- 
anism of nonperturbative gluon fluctuations which leads to confinement also generates an 
interaction of the strings in colliding hadrons. In this picture hadron-hadron scattering can- 
not be constructed from quark-quark scattering since the string-string interaction plays an 
important role. It leads to cross sections which are determined by the transverse extensions 
of the interacting hadrons. In addition to the forward scattering amplitudes, the model pro- 
vides the t-dependence of the cross section explaining the phenomenologically observed 
correlation between elastic slopes and total cross sections. 

Soft electroproduction on the nucleon can be calculated along the same lines as hadron- 
nucleon scattering using a model wave function for the photon. At small Q^, the photon 
is of hadronic size and large-distance-physics as in hadron-hadron scattering should apply. 
In another paper 0, we construct the wave function of the photon as a superposition of 
vector meson states and calculate the production cross section of p, p' and p" . This approach 
is limited to virtualities 2 GeV^ and the unknown couplings of the electromagnetic 

current to the p' and p" introduce new parameters. In this paper, we explore the possibility 
to represent the incoming photon as a qq state. With increasing Q^, the transverse extension 
of the qq dipole diminishes in a way that depends on the polarization of the virtual photon. 
We shall demonstrate how this mechanism shows up phenomenologically. This allows us 
to study the transition from large to short-distance-dominated-processes. The cross section 
for transversely polarized photons has a large non perturbative part, because endpoints at 
momentum fraction 2; = and 2; = 1 in the photon wave function do not select a qq system 
of small transverse separation. Therefore our large distance mechanism is important here, 
too. 

In our model the length of the string connecting the valence quarks in the hadron turns 
out to be very important. This length depends on the light cone wave function of the hadron. 
There is yet little knowledge about the physics determining the light cone Hamiltonian in 
nonperturbative QCD. So at the moment these wave functions and their integrated distri- 
bution amplitudes fulfill mainly a phenomenological task to parametrize the valence quark 
content of the hadron. Although, the exact value of cross section depends in our model on 
the detailed form of the wave function, the Q^-behavior of specific observables such as ratios 



2 



of longitudinal to transverse vector meson production or elastic slopes are likely to provide 
a good test for the string picture inherent in our model. 

Electroproduction of vector mesons has also been discussed within a soft pomeron frame- 
work in Ref. 0J^. In this model, transverse sizes of hadrons only play a marginal role be- 
cause, on the one hand, hadron scattering can systematically be reduced to quark scattering 
through the property of quark additivity of the model and, on the other hand, the vec- 
tor meson wave functions are assumed to be wider than the distance of the quarks in the 
virtual photon and are thus replaced by their value at the origin. These assumptions are 
phenomenologically tenable if one further assumes that the quark-pomeron coupling is fla- 
vor dependent. At intermediate Q"^ = 1 — 10 GeV^ 0, a pomeron form factor is used for 
far-off-shell quark legs. At larger values of [§, nonperturbative two-gluon exchange is 
applied which leads to color singlet cancellation at small gg-dipole size. 

In a series of papers 0, a perturbative two-gluon exchange model extended to include 
nonperturbative effects via the gluon distribution in the proton has been developped to 
evaluate vector meson production. This approach has been further refined to incorporate a 
BFKL-like evolution to accomodate both energy and Q^-dependence. Dipole scattering is 
the basis of this framework. In the following paper we are treating the photon and vector 
meson in a similar way. The main difference between our approaches lies in the reaction 
mechanism for soft diffraction. 

The importance of the gluon distribution as a necessary part of hard diffraction has 



been advocated in Ref. [p!0| , pT| . These authors have limited the range of applicability of 
their perturbative calculations to reactions where a large transverse momentum scale rules 
the exchanges, i.e. to heavy quark production such as J/'i/'-production |T0| or at large 
(Q2 > lOGeV^ fn]]). 



HERA has opened up new possibilities to enlarge the energy and range. Recent ZEUS 
and HI data indicate that at large the cross section may rise more steeply than expected 
from soft pomeron exchange. A possible explanation in the language used above [|T^] is that 
the evolution of the wavefunction at higher energy and gives rise to more and more 
dipoles inside the hadron. This evolution would then be a prerequisite to discuss the energy 
dependence of the virtual photon cross section with increasing resolution. This phenomenon 
of "hard" pomeron exchange will not be addressed in the following study which deals with 
soft pomeron physics at a fixed energy. 

The kinematics are defined in Fig. |I|. We denote the initial photon 4-momentum with 
q, the initial nucleon momentum with p, and the equivalent final states with q' and p'. 
A = q' — q is the momentum transfer and the independent Lorentz invariants are 

s = {p + g)^ 

t = A' = {q'~ g)^ 

We are interested in soft reactions, i.e. |t| < 1 GeV^, at high energy, s ^ and s ^ 
(e.g. s > 100 GeV^). In this domain, xb = Q"^ /2p.q is small, e.g. xb < 0.1. 

To be specific, we use the center-of-mass frame where the photon momentum points 
along the ^-axis. Then the absolute sizes of the 3-momenta are given as 
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p p 

FIG. 1. Kinematics of the reaction + p ^ V + p. 



2 2 s 



In this frame the vector meson emerges with a small transverse momentum |q(p| = At ~ 
6*^^/2 < 1 GeV. Unlike in the elastic scattering case, the center of mass momentum varies 
in this reaction 

X II I 'I Q' + ^F 

This implies that the momentum transfer A has a time component A° f=i 5, besides the 
space components ^ —S — ^s6'^/4 and At ^ \^ 9/2. Let us notice for completness that 
the square of A is t = A^ -A^ + to with to = -M^{Q'^ + M^f/s^. In the high energy 
hmit s > + the space components dominate. When we demand in addition a finite 
transfer, i.e. a small scattering angle 9/2 = 0{l/y/s), the transverse component is the leading 
component of A. In the following we shall therefore neglect all component besides At. 

The outline of the paper is as follows. Section II gives a very short and non-technical 

description of the model of the stochastic vacuum and describes its application to high-energy 
scattering. Section III deals with the specific features of electroproduction, i.e. the photon 
and vector meson wave functions, which are used in the evaluation of the longitudinal and 
transverse cross sections. Sections IV contains the numerical results for p, (p, J/ ip integrated 
and differential cross sections as functions of Q^. As far as possible these results are compared 
to experiment. Section V concludes with a discussion of the results. 



II. HIGH-ENERGY ELASTIC SCATTERING IN THE STOCHASTIC VACUUM 

A. The model of the stochastic vacuum 

The model of the stochastic vacuum is based on the asumption that the contributions of 
the slowly varying gluon fields in an infrared regular QCD can be approximated by a simple 



4 



stochastic process (for a review see Ref. [|T3|)- Already the assumption that this process has 
a converging cluster expansion leads to linear confinement in a non-Abelian gauge theory. 
As usual, approximations to a quantum field theory in the functional approach are more 
safely made in an Euclidean rather than Minkowskian field theory. It turns out however 
that in high-energy scattering there is no feasible way to continue Greens functions from 
the Euclidean to the Minkowskian world and we therefore have to formulate the model in 
the Minkowski continuum. This seems at first sight more dramatic than it turns out to 
be finally, since at the end we have to evaluate the relevant quantities only at spacelike 
Euclidean distances, i.e we can take these quantities from an Euclidean field theory. 

In order to define gauge invariant correlators we introduce the modified gluon field 
strength F^y{x, uo) which is obtained from the field strength at point x by parallel transport- 
ing the colour content to the point u 

Ff,^{x,uj) = ^~\x,uj)Ff,^{x)^x,uj), 

with cj) = Pexp[—ig Adz]. 

Assuming that the main features of the correlator {F^^{x,uj)Fp„{y,u!)) do not depend 
crucially on the choice of the reference point u we obtain for the dependence on z = x — y 
the most general form 



{g^F;^{x,uj)Ff^{y,u))A = ^^^-(/FF){«:(r/^,r/,. - r/^.r/.p)D(zV«') (2-1) 

+ (1 - K)]^[df,{zpr]^^ - z^r]^p) + d^{z„r]pp - Zpfj^^) Di{z^/a^)Y 

The correlator D is typical for a non-Abelian gauge theory (or an Abelian theory with 
monopoles) since the homogeneous Maxwell equations 

e^^'^'d^Fp^ = 

allow only the tensor structure proportinal to Di, hence k = in an Abelian theory without 
monopoles. 

In a Gaussian model, where all higher cumulants in the linked cluster expansion fl^ 



are 



neglected, we obtain a relation between the slope of the static quark-antiquark potential and 
the typically non-Abelian correlator D 



^ /■ + 00 

a = K—(g^FF)a'^ duD(-u). 
144^ 'Jo 



The choice of phenomenological parameters will be given in Sec. II C 



B. Scattering of two color singlet dipoles 

The high-energy scattering of two color singlet dipoles gigi and ^2^2 can be treated anal- 
ogously to the situation of heavy quarks encountered in the Wilson area law. The relativistic 
quarks and antiquarks move along two opposite straight line trajectories on the light cone. 
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FIG. 2. Configuration of the two interacting loops in the transverse plane. With our choice of 
frame, the loops 1 and 2 lie, in the (x'', x^)-plane, on the lines = and x^ = — respectively. 



In order to apply the model of the stochastic vacuum to high-energy hadron hadron scat- 
tering we adopt the method of Ref. 0. In this approach the problem is first considered as 
the scattering of quarks in an external colour field which is solved for fast moving quarks by 
the leading term of an eikonal expansion, i.e. the quark picks up the eikonal phase 



V = exp[—ig J Adz] 



where F is the classical path of the quark. 

This phase is manifestly gauge dependent, but if we consider a fast moving dipole, i.e. a 
quark and an antiquark moving on parallel lightlike trajectories connected by a Schwinger 
string, then we have to evaluate rather a Wilson loop 

W = exp\-ig / Adz], 

JdS 

than the path integral above. The open ends of the q and q trajectories in dipole 1 and dipole 
2 are closed by small transverse lines yielding two loops, dSi and 882, which have transverse 
extensions according to the lengths of the dipoles ri and r2. The dipoles are positioned 
relative to each other with a given impact parameter b. The loop-loop interaction amplitude 
for this system of dipoles is calculated in Ref. M 



J(xi,Xl,X2,X2) 



^_Ltr (iyi(xi, xi) - 1) ^tr (iy2(x2, X2) - 1) 



where the bold faced vectors Xj = {xl,xj) denote the two-dimensional positions of quarks 
z = 1, 2 in the transverse plane {% refers to the corresponding antiquarks). The geometry of 
the loops is shown in Fig. 0. 

The various steps and approximations necessary to derive a tractable expression of J 
have been developped in Ref. P| to which we direct the interested reader. Here we only 
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summarize these steps. First one transforms the hne integrals appearing in the non abehan 
phases Wi and W2 into surface integrals. The manipulation of gauge invariant quantities 
leads to the introduction of a reference point u in between the two surfaces. The surfaces 
to be considered are two pyramids, 5*1 and S2, with u as apex and dSi, 882 as respective 
basis. This is shown in Fig. ^. The surfaces are the world sheets of infinitely many gluons 
in the two-hadron state which interact via the correlator Eq. In order to make the 

calculation practical, these interactions are truncated to fourth order in the field strengths. 
Because of the Gaussian process the terms factorize into products of integrals over the 
correlators linking the surfaces generated by Wi and W2- One gets 

Symbolically, this expression can be rearranged into a sum of four interaction-terms, x? 

" 8Nl;{Nlj - 1) ( 1^^^ ) + x{qiq2) - x{qiQ2) - x(glg2)}^ (2.2) 

where each term has confining and non-confining parts from the basic gluon-gluon correlator 
in the vacuum, xill) = i^Xcill) + (1 ~ i^)Xnc{'i<l) , which we now specify. 

It is possible to choose functions D and Di in a way that the natural fall ofF at large 
Euclidean distances does not impede a meaningful continuation to Minkowski space. A 
possible choice consistent with confinement and lattice computation is 



The function Di is a priori an independent one, however lattice simulation indicates that 
the behavior of D and Di is similar and we shall assume for simplicity that Di = D. 

For the non-confining part the two surface integrals can be performed and result in 
quark-quark interaction terms such as 




Xnc{qiq2) = 

Since nonperturbative gluon correlations are of size a, two color charges can only interact 
if their trajectories enter in a common domain of size a. A constituent quark picture arises 
where the elementary color charges are surrounded with gluon clouds. This can be seen 
in Fig. Ka) where the interaction amplitude J for k = between a dipole target of size 
r2 = 12 a oriented along a given x-axis and a dipole probe of size ri = a is plotted as a 
function of the impact position. For simplicity we sum over the orientation of the probe. It 
turns out that with a physical correlation length around 0.3 fm a physical target has a size 
of 4 a-5 a so that constituents in the target are not as well separated as in Fig. ^ 

For the confining part on the contrary the integrals have a path dependence which is 
linked to the non-abelian nature of the confining term. Physically this means that the color 
dipoles connected by their strings interact as whole objects rather than as isolated endpoints. 
In space-time the integration over the surfaces is done with the reference point u chosen in 
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FIG. 3. Space-time representation of the pyramids and their transverse projection. The shding 
sides of the pyramids give the domains and ^2 of the surface integrations. Note that the two 
Wilson loops are not parallel in the transverse plane. 
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FIG. 4. (a) Color interaction amplitude Eq. ( |2.2| ) for the non-confining case, k = 0, as a function 
of the impact position between the two dipoles. Dipole 1 has a transverse size ri = a and we sum 
over its orientation. Dipole 2 has a transverse size r2 = 12 a and lies along the x-axis. (b) Color 
interaction amplitude for the confining case, k = 1. 
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the most symmetrical way as shown in Fig. |^. The result depends only very weakly on the 
choice of the reference point and has the following form: 



(2.3) 

K2 (^^|xi -x^ -a(x2 -x<,)|j + (1 ^ 2)| . 

The angle ipu denotes the angle between the vectors xi — x^ and X2 — x^^. The amplitude in 
the string-string interaction picture, i.e. k = 1, is shown in Fig. |^(b) with the same choice 
for the target and probe sizes as in the non-confining case. The interaction is non-zero 
whenever the probe is close to the line connecting the target quark and antiquark. We note 
that the string and constituent picture both differ from the optical droplet picture where 
the charge distribution form factor is responsible for the differential cross section. Via the 
wave functions of the valence quarks the geometrical sizes of the hadrons enter in the cross 
sections. 



TT 



Xc{qiq2) = 2 cosV^i2 



da 



Xl 



a(X2 



a 



C. From dipole-dipole to hadron-hadron cross section 

A valence quark picture can be constructed from nonperturbative scattering amplitude of 
color dipoles with fixed lengths ri and r2 by distributing the positions of the end-points of the 
strings according to a quantum mechanical wave function. Since for high-energy scattering 
the incoming particles propagate along the light cone, it is natural to choose light cone wave 
functions. The amplitude of the process can be written as 

M = 2^sJ ^^^U,{z„r,) J ^^|^,(.2, r2)| V({x.}), (2.4) 

where the index "1" refers to the photon or vector meson side whereas the index "2" is 
attached to the nucleon coordinates. (Conventions are fixed in Appendix ^.) For simplicity, 
the nucleon is considered in a quark-diquark configuration. It has been shown in Ref. [0] that 
quark-diquark and three-quark pictures lead to similar predictions for diffractive scattering 
once the model parameters are adjusted to fit the proton-proton cross section. In the fol- 
lowing, we also fix the parameters to fit proton-proton cross section and therefore we do not 
expect any significative dependence on the model of the proton. For C = P = —1 exchange, 
the quark-diquark is favoured since it suppresses the odderon contribution |r7|. 



The ip^s are the valence light cone wave functions of the corresponding hadrons. They are 
usually defined in momentum space where they describe the probability amplitudes to find in 
a hadron with momentum {P~^, P}Q and well defined angular momentum and flavor content 
a quark and an antiquark with momenta {zP^, k+zP} and {{l — z)P'^, — k+ (1 — 2;)P}. One 
crucial property of light cone wave functions is their dependence on z and k alone |jl8|,|l9|- 



^ Light cone coordinates are = {P^ it P'^)/\/2, P = = {P^ ^ P^). A particle momentum is 
fully specified by the set {P+, P} (P" = (P^ + M2)/2P+ with M the particle mass). 
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This implies that upon Fourier transformation in the transverse plane, the relative coordinate 
of the qq pair, r = — x^, is easily separated from the position of the hadron "center" , 
X = z^q + (1 — z):x.q. In a similar way to non relativistic physics, the non trivial degree of 
freedom of the wave function can be isolated, the resulting transition matrix element is 



iP X 



(z,Xg,X5|P+,P) = ijj{z,r)e 

In the amplitude given in Eq. ( |2.4| ) the impact parameter b denotes the transverse 
separation between hadron centers, b = Xi — X2. The transverse positions of quarks are 
then given by 

xi = xo + b/2 + (1 - zi)vi, 
X2 = xo - b/2 + (1 - Z2)r2, 
xo= (Xi + X2)/2, 

i.e. xq is in the center of Xi and X2. Antiquark positions are related to quark positions as 



Xj = Xj — Fj. We notice that the reference point uj chosen in Sec. [II B| "moves" with respect 
to Xq, Xj^ = Xq — l/2[(zi — l/2)ri — (^2 " l/2)r2]. We have, however, checked that choosing 
Xo rather than uj as the reference point for the computation of the loop-loop amplitude has 
negligible numerical effects. 

In the present study, we are interested in electroproduction of different vector mesons 
under various kinematical conditions on a fixed proton target. It is therefore instructive 
to isolate the variable part of the amplitude associated with the photon and vector meson 
coordinates, -Zi, ri, by integrating out the nucleon coordinates, 2:2, i'2- To this end we define 

Jp(2i,ri, At) = 2j^ bdb27iJo{ATb) J -^^\M Z2, 1^2)1 J{b, zi,ri, Z2,r2), (2.5) 
so that the amplitude Eq. (|2.4|) is now written as 



/dz r 
— ^^—^'ipli)^{zi,vi)Jp{zi,vi,AT). (2.6) 

The determination of parameters of the model of the stochastic vacuum can be made 
in different ways. We follow the method given in Ref. We use as input parameters 
the total proton-proton cross section at ^/s = 20 GeV, namely app = 35 mb and the slope 
Bpp = 11.5 GeV~^ of the p-p elastic cross section. From lattice simulations [^, we take the 



mixing coefficient k = 0.74 and the curve relating gluon condensate and correlation length. 
The square of the proton quark-diquark wave function is taken in the simple form 

\^{z,r)\^ = 4uj^6{z - l/2)e~^"'\ 

As ouput we obtain the gluon condensate (g'^FF) = 2.49 GeV^, the correlation length a = 
0.346 fm and the proton transverse radius Rp = 1.51a = 0.52 fm. The parameters are different 
from those of Ref. where the influence of the non-confining term, Di, in the correlator 
Eq. ( p.l|) was neglected, i.e. Di was set to 0. The proton size is smaller than the transverse 
rms charge radius, Rt = 0.68 fm. Sea-quark contributions may enter into the form factor 
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increasing the charge radius compared to the radius of the valence quarks. The same sea 
quarks or extra color dipoles may also enter in the cross section at higher energies where 
it is increasing with energy. New lattice simulations try to isolate the perturbative 
contributions from the nonperturbative gluon fluctuations in a more precise way and may 
be incorporated together with a better treatment of the perturbative two-gluon exchange. 




1 2 3 4 5 r^{fm^) 

FIG. 5. Dipole-proton total cross section as a function of the dipole size for z = and z = 0.5. 
The dependence in z is rather marginal as it becomes noticeable only for very large separation 
of the qq pair. The cross section behaves as aqq oc r" with n = 2 for small extension and slowly 
decreasing at larger distances. 

In Fig. ^, we show the behavior of the function 

/■27r lift 

4'\z,r,AT = 0)= —J,{z,r,e,AT = 0). 

It represents the total cross section of a qq dipole of fixed size r averaged over its orientation. 
For varying dipole sizes the total cross section on the proton increases quadratically for small 
dipoles until a size of r = 1-2 a, then the increase continues but with a decreasing power. 



This feature is distinct to model of perturbative gluon exchange in Ref. |^ where this total 
cross section saturate at about twice the proton radius. 

It is also instructive to consider the A-r-fall off of the transition amplitude on the hadron 
radii and the correlation length a. In Ref. [Q, the logarithmic slope B of elastic cross sections 
has been numerically parametrized as 



B ^ 1.56 0^ + 0.24 {Rl + Bj. 



The large first term which is independent of the radii is specific for the model. 
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III. ELECTROPRODUCTION OF VECTOR MESONS 



In this section, we describe the qq wave functions of the photon and vector meson which 
enter into the expression of the production amphtude Eq. (^I6|) . Even at the quahtative level 
there are still large uncertainties about the correct dynamical description of hadrons. Indeed, 
one of the most interesting issues of current and future experiments is to shed light on the 
long distance properties of QCD. This evidently includes the unravelling of basic facts about 
hadronic wave functions. We shall show that the electroproduction of vector mesons is quite 
sensitive to their wave function. 

For the time being, one has to make assumptions which influence the result as strongly 
as the dynamical features of the transition operator for diffraction given by the model of 
the stochastic vacuum. In the analysis of hadron-hadron scattering a simple transverse 
wave function for the vr meson has been chosen 



Note that in this reference the dependence of the geometry of the loop-loop interaction 
on the respective light cone fractions of the quarks and antiquarks was not yet considered. 
Including this dependence and wave functions ip{z, r) with a reasonable behavior on the light 
cone momentum fraction z, we found similar results to those given in Ref. 0. The transverse 
size of the studied hadrons still determines the size of hadron-hadron cross sections. We shall 
see in the following how this feature is modified when one is considering photon induced 
reaction. 



The qq wave function of the photon carries as labels the virtuality and the polarization 
state A of the photon. It describes the probability amplitude to find a quark- ant iquark 
pair inside the photon with light cone fractions {z, 1 — z) and transverse separation r = 
(r cos ^, r sin ^). The qq state is in a configuration with given flavor (/, /) and helicities 
{h,h). The color part of the wave function is treated separately and considered together 
with the Wilson loop in the way described in Sec. [II B| and we are only left here with an 
overall factor y/N^. The photon couples to the electric charge of the quarks with ef6fj where 
6/ = 2/3 e or — l/3e respectively. The helicity and spatial configuration part of the wave 
function looks different for various photon polarizations. It can be computed in light cone 
perturbation theory and one has to lowest order (see the Appendix ^ 



(3.1) 



A. Photon wave function 




(3.2) 



with 





V'7(Q2,0) 




V2 [ie'e[z5h+8j,_ - (1 - z)5H.5n+) 
V2 Ue-''e({l - z)5h+h_ - z6h-h+ 




Kijer) 
2tt 
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where e = yz{l — z)Q'^ + w?j and nif the current quark masses given in Table | for the 
different flavors. Kq, Ki are modified Bessel functions. 

The longitudinal photon wave function is peaked around z = 1/2, so that the longi- 
tudinal photon interacts like a small dipole, r ~ 1/Q, at large Q^. On the contrary, the 
transverse photon is almost fiat in z so that, at large Q^, it interacts partly like a small 
object for intermediate z and partly like a large one for z ~ ruf/Q when light quarks are 
involved. For heavy quarks, c, b, the inverse of the quark mass limits the photon extension. 
In electroproduction of vector mesons the effective dipole size is fixed by the overlap of wave 
functions of the photon with the vector meson. Due to the shape of the latter the small z 
region is somewhat suppressed and the transverse region explored for Q above 1-2 GeV is 
below or around 1 fm. Lacking a better knowledge in the region of large transverse size, we 
use the above wave function. This may be tested by forthcoming experiments on the ratio 
of longitudinal to transverse cross sections. 

B. Vector meson wave function 

For the hadron wave function, we use the same notation as described above for the 
photon. We have already taken care of color in the construction of the Wilson loops cf. 
Sec. |II B| . For the flavor content we consider that the p°(770), ct;(782), 0(1020) and J/ip 
mesons are respectively pure isospin 1, isospin 0, ss and cc vector mesons. This flavor part 
together with the configuration part ip forms the wave function ipv needed in Eq. (|2.6| ). 

Some information has been obtained on the spatial dependence of hadron wave functions. 
The first piece of information comes from long distance physics which various quark mod- 
els describe successfully. These models tell us that a hadron at rest can be modelled with 
Gaussian wave functions as a sytem of constituent quarks moving in a harmonic oscillator 
potential. To boost non relativistic wave functions to a fast moving system is not a trivial 
step. Technically the interplay between the transverse and longitudinal dynamics in light 
cone physics as well as the treatment of spin degrees of freedom remain to be understood. 
Within the model of the stochastic vacuum a nice quantitative description of hadron-hadron 
soft collisions has been obtained by disregarding spin and light cone fraction dependences. 
Therefore we assume that for soft collisions between large objects, a simple-minded descrip- 
tion of hadrons Eq. ( |3.1| ) suffices. A smooth 2;-dependence on the light cone momentum 
fraction will not change this picture. 

The second piece of information comes from short distance physics and perturbative 
QCD supplemented by sum-rules, where some properties of the valence wave function at 
0-transverse separation are known (in particular the end-point z — *■ 0, 1 behavior can be 
analyzed). When wave functions at short distances are involved, e.g. in hard exclusive scat- 
terings |]22|, the value of the meson wave function at the origin determines the value of cross 
sections. The wave function at the origin is related to the leptonic decay width of the meson. 

The above observations allow to make an ansatz for vector meson wave functions: 

ipvio) = - z)^j^ f[z)e 2 , 
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This ansatz has the following properties. The main transverse dependence ^-'^^'''^1'^ and the 




function f{z) are modeled in the way proposed by Wirbel, Stech and Bauer 



f{z) = M^z{l-z))e-'''vi--^/^f/^-\ (3.4) 

We use the same functional form for all vector mesons. The transverse size parameter uj 
is related to the vector meso n rad ius, = {x^ + y"^ + z^)-*^/^, which for the Gaussian 

shape adopted is {R^Y^"^ = \j3/2/2u. This quantity is presumably not very different from 
the electromagnetic radius which unfortunatly is unknown for vector mesons. The way out 
in the quark model is to fix u and J\f by the normalization and the e~^e~ decay width (see 
Appendix 0). We draw attention to the fact that applying this procedure to the above 
parametrization of the wave function leads to different sets of parameters, {uj,Af}, for longi- 
tudinal and transversal mesons. We find radii in this way which are reasonable in the whole 



family of vector mesons. This is a welcome property because, as we shall see in Sec. the 
vector meson transverse size is one of the important ingredients determining the cross sec- 
tions in the intermediate Q^-range. The above form is written to have explicitly the correct 
value of the wave function at the origin, i.e. the gg-state fulfills the equation 

{0\riO)\Viq,X)) = efvMvs'^iqA), 



with fv the meson decay constant. In Eq. ( |3.3| ), iy is the mean quark charge in the meson 
state in units of the proton charge (see Table H in Appendix IB]). The helicity dependence of 
the wave functions is modeled after the perturbative 7 — > gg transition. 



C. Cross sections 



Let us collect the formula necessary to compute cross sections. It is convenient to expand 
the quantity Jp in Eq.(^.5D in terms of L^-eigenfunctions 

Jp{z,, n, e,, A) = ^ e^-^Vf )(;2i, ri. A), 

m 

where thanks to the periodicity of Jp only even m are present. It follows that in the reac- 
tion studied, only transverse to transverse and longitudinal to longitudinal transitions are 
expected. From the explicit form of the wave function, Eqs. (|3.2| ) and ( |3.3| ), there is a pos- 
sibility of helicity change by two units in the process. We have, however, observed that the 
corresponding contribution to the cross section is smaller than 2% in the whole Q^-range. 
We disregard this contribution. 

In the following, we distinguish transverse and longitudinal cross sections which, in the 
present conventions, are 
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daL _ 


1 


dt 


167r 


dax 


1 


dt 


IGtt 



dziTidvi X 
2 ^^(0) 



dziVidri X 

in I 



(3.5) 



where the average over proton hehcities is understood. The combination '?/'|/(a)V^7(a) (-^i, ''^i) are 
computed by multiplying Eq. ( |3.2| ) by Eq. ( ^.3[ ) supplemented by the flavor part described 
in Sec. [Ill B| giving 



u eri r 



ml 



(3.6) 



Experimentally, the differential cross sections Eq. ( p.5| ) are difficult to measure. On the 
one hand, the separation of transverse and longitudinal cross section is not easily done, on 
the other hand accurate data exist only for t-integrated cross section. Nevertheless, some 
results have been obtained for 



da, 



exp 



daL dax 



dt dt dt 

as a function of A^. e is the rate of longitudinally polarized photons which depends on the 
lepton scattering angle, 6'e, and the photon energy, u. In the proton rest frame 



1 + 2 (l + z/VQ^) tan2(ee/2) 



-1 



We shall also compare our theoretical results with the integrated cross section cTexp = cctl + 
(Tt, for various Q^. By analyzing the vector meson decay, it is possible to check the validity 
of s-channel helicity conservation or assuming helicity conservation to deduce R = aL/aT- 



IV. NUMERICAL INVESTIGATION 
A. General results 

Before entering in a detailed analysis, let us discuss qualitatively the Q^-dependence 
of differential cross section expected from Eq. (|3.5|). The quantity J^\zi,ri,t) is power 
behaved in ri and decreases exponentially with \t\. It depends only weakly on zi. At fixed t, 
J^'^^ oc rf with a = 2 at small ri and a slowly decreasing for ri ^ a (see Fig. Given this 
behavior, we also need an estimate of the effective size of the photon-vector meson overlap 
after integration over zi in Eq. (|3.6| ). One can anticipate two extreme regimes. At small Q^, 
the effective size is driven by the meson wave function so that daL oc and dax is constant. 
At large Q^, the effective size is given by the photon wave function alone, which leads to 
daL oc Q^^ and dax oc Q^^. 

These asymptotic behaviors of cross section are rather model-independent. Most of the 
current experimental results, however, are below the asymptotic large Q^-region. The the- 
oretical behavior of cross sections in our calculation shows a much more specific depen- 
dence on in this region. In Fig. ^ we show the logarithmic p-production cross section 
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da/dt{t = 0) (nb.GeV-2) 
10000 F 1 




1 10 g'(GeV2) 

FIG. 6. da/dt{t = 0) for p-production as a function of for the fongitudinal (full) and 
transverse (dots) cross section. The effective power of the fall off with Q is increased by 2 units in 
the range 1-10 GeV'^ but the asymptotic behavior is only reached in the 10-100 GeV^ 

\ogda / dt{t = 0) as a function of Q^. The effective power n of the fall off da /dt{t = 0) oc 
is about 2.5 and 4 at = 1 GeV^ for the longitudinal and transversal cross section respec- 
tively and this power increases to about 4.5 and 6 respectively at = lOGeV^. Thus the 
asymptotic regime starts only above that value. 

As we have emphasized the result depends on the wave function chosen to represent the 
vector meson state. Without changing the parametrization, this can be seen by varying the 
size parameter uj of the vector meson wave function. With the value of the wave function 
at the origin kept fixed we observe that decreasing uj by 5% increases the cross section at 
Q2 = 1 GeV^ by about 20%. This modification becomes less than 3% at = ig GeV^ 
because the vector meson size is less important for a small q-q state in the photon. One can 
also think of changing the light cone distribution of the quark and antiquark in the meson 
wave function. An indication of the resulting modification is given by changing the factor 
\l z{l — z) — > [z{l — 2;)]^/^ in the parametrization of f{z) in Eq. ( |3.4| ). This leads to about 
30% decrease of the cross section in the whole range examined. 

In Fig. 1^, we show the importance of the transversal extension of the virtual photon. We 
introduce a transversal cut-off in the cross section, i.e. we cut the ri-integration in Eq. ( p.6|) 
at a fixed value ri < rj"*. As expected the gg-wavefunction in a transversal photon extends 
to much larger values of ri, e.g. at = 4GeV^ more than 50% of the cross section comes 
from transverse separations bigger than 1 fm. In the longitudinal case, the contribution of 
the region with transverse separations ri > 1 fm drops down to 15% at = 4 GeV^. 
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0.5 1 1.5 2 2.5 r™* (fm) 




0.5 1 1.5 2 2.5 rj"* (fm) 

FIG. 7. (a) The contribution to the transverse differential cross section of p-production with 
respect to the radial parameter of the photon and vector meson side. The integration over the 

transverse distance ri is taken up to a cut-off r™* for several values of the photon virtuality, 
= 1, A and lOGeV^. The transverse differential cross section at t = is plotted as a function 
of rf"*. It is normalized to its value with the cut-off removed, (b) Same as (a) for longitudinal 
polarization. With increasing the cross section is dominated by short transverse distances in 
both cases but the saturation occurs earlier for the longitudinal cross section. 
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FIG. 8. The scaled cross section (t{Q'^) for p-production in nb.GeV^. The circles are the 
NMC-results [^] and the diamonds represent our prediction for the quantity Q'^ (ccl + cj-) with 
the experimental polarization rate of NMC. 



B. p, Lo, (^-production 

In the range from 1 to lOGeV^, the behaviour of the production cross section 

observed by EMC and NMC is very well reproduced by our calculations. Besides the Q^- 
dependence, also the absolute values of the cross section are reproduced. It should be noted 
that within our model we have introduced no new parameters and the parameters under- 
lying the interaction on the quark-gluon level are determined by soft high-energy proton- 
(anti)proton scattering. We show in Fig. || NMC deuteron data together with our pre- 
diction for (5'*(o"l + (Jt). We notice that this theoretical quantity is not the actually measured 
cross section due to the polarization rate of the photon, e 7^ 1. Therefore we also show with 
diamonds the quantity Q^{eaL + (Tt) at the Q^-points of NMC using their value of e(Q^). 

In our approach, the approximate behavior is due to a combination of different fall 
offs in (Tl and ax which themselves come from the interplay between the size dependence of 
the dipole-proton cross section and the effective size of the photon-meson overlap. This is 
very different from the dynamics which occur in quark-quark scattering which leads however 
to a similar Q^-dependence. As explained above the asymptotic behaviour for large is 
just governed by the dipole size of the virtual photon and thus model independent. 

A possible way to distinguish both approaches would be a precise measurement of the 
Q^-behavior of the ratio of longitudinal to transverse cross sections RiQ"^) = ctl/<^t- We 
plot this ratio in Fig. ^. For large Q^, R oc but this behavior is not yet reached in the 
intermediate range where R grows slower than Q^. Here we expect a rather different behavior 
in our model compared to other models. 

Another important check is provided by looking at the pr-dependence of the differential 
cross section. We show our result for edai/dt + dax/dt versus at 6GeV^ and compare in 
Fig. ITO| to the NMC points for the deuteron outside of the coherent production region . 
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2 4 6 8 10 Q'(GeV2) 

FIG. 9. The ratio for longitudinal to transverse cross section for p-production. The data point is 



from Ref. |24]. Other data compare well within errors but they either are far outside the 10-20 GeV 
range or have large err or bars. 



Notice that the t-dependence of our computation is not fully exponential. 

There is a non trivial Q^-dependence of the slope B in our model which is due to a 
decreasing transverse region probed by the slowly shrinking size of the photon as grows. 
We show this in Fig. |r^. 

For the production of the cu-meson, we expect a wave function very similar to the one of 
the p and correspondingly the ratio should be determined through the flavour factor 

fl/fl ^ 9%. 

Indeed, this is observed. 

The situation is more complex for the production of heavier vector mesons where the 
different quark content has several consequences. The direct appearance of a mass term 
in the meson wave function gives an additionnal component in the overlap Eq. (|3.6|) and 
also modifies the photon extension parameter = z{l — z)Q'^ + m'j. The quark content 
also influences both the transversal size of the vector meson and its momentum fraction 
distribution f{z). Whereas the first effects are quite easily controlable, nothing precise is 
known about the quantitative effect of a heavier mass on the distribution f{z). Qualitatively, 
one expects that the distribution becomes more peaked at 2; = 1/2 as the mass of the 
constituents increases. 

For 0-meson, a smaller transversal extension than in the p-meson is expected due to 
the heavier s-quark mass. The smaller size of the 0-meson reduces the 0-production cross 
section at low Q'^ values beyond the flavour factor 

fl/f^ ^ 27%. 

Such an effect is also observed in the difference between pion-nucleon and kaon-nucleon 
scattering in the model M . It should diminish with increasing in electroproduction since 
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0.1 0.2 0.3 0.4 0.5 0.6 Af, (GeV^) 

FIG. 10. The differential cross section, da/dt{A^), for j* + p ^ p + p at = GGeV^. Data 
from Ref. [||. 
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FIG. 11. The differential cross section for p-production, da/dt{Aj^), for a longitudinal photon 
(full lines) and a transverse one (dotted lines) at = 2 and lOGeV^. 
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FIG. 12. i?i>-production cross section compared to NMC data \24\. In our model the difference 
between theoretical and experimental result is attributed to wave function effect. 



then the amplitude is less and less determined by the extension of the produced meson, but 



rather by the virtual photon. This effect may have been observed in the ZEUS data . Also 
the change of the longitudinal distribution in z between and p-meson may influence the 
cross section independently of Q^. This difference is absent with our choice of distribution 
for the 0-meson which is numerically the same as for the p-meson. Our resulting theoretical 
cross section for the 0-production reproduce the dependence of the NMC data, but its 
absolute value is practically a factor 2 too large (see Fig. [12|) . 

The increase of the longitudinal to transverse ratio, R{Q'^), looks the same as for the 
p-meson, its overall magnitude being just reduced by about 20%. The pr-dependence also 
exhibits a similar pattern with a small broadening of the diffraction peak as is increased. 



C. J/V'-production 

For heavier quark pairs, the large quark mass leads to more dramatic modifications. Let 
us first notice that in our model we assume that the quarks move on lightlike trajectories. 
This can only be the case at energies far above 2 mj. Therefore a center-of-mass energy bigger 
than 10 GeV is necessary in the J/ip-case. At these energies a moderate energy dependence is 
observed which is not contained in our model. In return, the large quark mass provides a hard 
scale so that also photoproduction data are accessible within our perturbative treatment of 
the photon. Strictly speaking, the difference to = ~Mp{Q'^ + MyY / s^, between t and — A^, 
leads to a phase space threshold, e^^° , for the J/-?/'-production in the energy range we are 
considering. At ^/s = 15 GeV, it is easy to see that this threshold effect is only sizeable for 
large e.g. for = 10 GeV^ and S = 5 - 10 GeV'^ one gets e^*" = 0.97 - 0.94. We shall 
disregard this factor in the following. 



The discussion of the asymptotic regime given in Sec. [IV A| can be refined in the presence 



of a large quark mass. Let us reexamine the short distance regime in the presence of the mass- 
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[l + QV^j]'^^(nb.GeV4) 

"1 1 1 r 




2 4 6 8 10 Q^GeV) 

FIG. 13. J/V'-production cross section for longitudinal (dots) and transverse (dashes) po- 



larizations. To compare with EMC-data [26|, one has to combine these two cross sections into 
a = eaL + (JT with the polarization rate measured by EMC, e 0.7. 



terms in Eq. (|3.6| ). To simplify further, we temporarily assume a simple non-relativistic form 
for the distribution, i.e. f{z) oc S{z — 1/2) and for consistency take = Mj/2. We consider 



the domain of large enough 
exp(— u;^r 2/2) ^ 1 and J(°)(;zi = 



= ml + QV4 = 
1/2, n, A = 0) 



dt 

dax 
IT 



Op 



(t = 0) 
(t = 0) = a. 



m 

« Crj 



J + Q )/4, where we can approximate 



One gets 



(MJ + 

(mJTW 



'MJ + 



From these expressions one sees that the relevant scale for J/'?/'-production is Mj + 
rather than Q^. As in the light quark case, the longitudinal cross section is expected to 



dominate at large Q^, namely the quantity (Mj + 8uj'^y/{Mj + 2uj^)^ ^ 1.4 leads to a 
ratio R ~ 0.7 Q"^ /Mj. The differential cross section is expected to fall off as da/dt oc (1 + 
eR){Mj + Q2)~^. Experimentally a Q^-dependence such as {Mj + Q^)""" is observed with n 
around 2 but the accuracy of the data is not sufficient to exclude a complete short distance 
fall off. Let us stress that an accurate test of this power-law is a necessary prerequisite to 
understand the physics at work. 

In Fig. 0, we plot (1 + Q'^/Mj)^ a, i.e. the cross section rescaled by the observed data 
fall off, for the transversal and longitudinal cross section separately together with the data 
recorded by EMC in the energy range y/s =10-20 GeV ||2^. The photoproduction comes out 
fairly in the energy range ^/s =10-20 GeV, where several measurements have been performed 
leading to a production cross section between 10 and 20 nb. As could be guessed from the 
study of the asymptotic behavior above, the shape of the 2;-distribution and the value of 
the charm quark mass determine the size of the cross section. Changing the quark mass 
by 5% would lead to a 20% change in the cross section at = and to a 10% change at 
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FIG. 14. da/dt{A'^) for J/V'-production at = q and lOGeV^. The upper curve is our 
prediction for the photoproduction differential cross section. It can be compared to the measurement 



of Ref. 1 27 1 (filled circles) and to the extrapolation to = of the EMC data |26|. Similar 
data have been measured by NMC |28| at = 1.5 GeV^. Also shown are the differential cross 
sections at = 10 GeV^ for longitudinal (dotted line) and transverse (lower full line) photons. 
The dependence of the slope on and polarization is marginal in the J/V'-case. 



= 10 GeV^ respectively. As increases, our expectation follows qualitatively the pattern 
depicted for the short distance regime although quantitatively the intermediate transverse 
distance somewhat contributes to give a fall off flatter than the short distance one. 

We next turn to the study of the py-dependence shown in Fig. 0. We find a good 
agreement with the photoproduction measurement and with the extrapolation of EMC 
(open circles). We also note that, contrarily to the "large" hadron case, there is practically 
no Q^-dependence of the pr-fall off in the J/ip-case. 



V. CONCLUSION 

We have calculated the longitudinal and transversal differential cross sections for diffrac- 
tive production of p, u, cf) and J/^/'-mesons in the range of 2 GeV^ < < lOGeV^. The 
hadronic part of our calculation is based on a model for nonperturbative QCD, the model 
of the stochastic vacuum |jl| . The parameters of the model which gives a unified description 
of low energy and soft high-energy scattering phenomena can be obtained from a variety 
of sources: hadron spectroscopy, QCD sum rules, high-energy proton-proton scattering and 
lattice calculations of the fundamental gluon field correlator. We have used a consistent pa- 
rameter set very similar to the one used for hadron-hadron scattering in Ref. 0]. The virtual 
photon and hadron wave functions are light cone wave functions motivated by perturbation 
theory for the photon and relativistic quark models for the hadrons, we thus do not have 
adjustable parameters. A specific feature of electroproduction is the dependence of the cross 
sections on the photon virtuality which is reproduced by the model almost perfectly. 
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Even at = 10 GeV^ the cross sections have not yet reached their asymptotic be- 
haviour. Our calculation is consistent with the observed ratio of longitudinal to transverse 
cross sections. Here precise data at different values of could discriminate between differ- 
ent models. In our model, we also can calculate the dependence on the tranverse momentum 
transfer and it nicely reproduces the available data. 

Depending on photon polarization, the range 2-10 GeV^ corresponds to an effective 
transverse qq size lying between 0.5 fm and 1.2 fm. This region just interpolates between the 
short distance domain and normal hadron diameters, thus allowing a natural extension of the 
phenomenology of hadron-hadron scattering where the model has been applied originally. 
Model dependent features of the light cone wave functions enter into the magnitude of the 
production cross section. A good experimental separation of ai and ctt can help to obtain 
a real breakthrough in our understanding of diffractive electroproduction, since the physics 
of the qq pair state is so much different for both photon polarizations. The extension to 
photoproduction and low electroproduction necessitates a modification of the simple 
perturbative qq wavefunction in the photon. 

The model of the stochastic vacuum cannot predict the behaviour of the cross sections 
as a function of the cm energy ^/s. If all parameters are fixed, it yields constant cross 
sections. We therefore have confined ourselves in this paper to the energy range 10 GeV< 
a/s < 20 GeV. The absolute values for the cross sections are correctly reproduced for p, 
uj and J/^/'-production. Our results for 0-production are about a factor 2 larger than the 
NMC data, we can, however, explain the Q^-dependence of the ratio of to p production 



as observed by ZEUS ||25 



In the present model the energy dependence of the total hadronic cross section and 
the slope of the elastic cross section can be obtained consistently by increasing the hadron 
radii (slightly) with energy. Due to the expected difference of the energy dependence of the 
perturbative and nonperturbative contributions (hard and soft pomeron) one has also to 
take into account perturbative contributions if one wants to obtain a realistic model of the 
energy dependence of diffractive electroproduction. This will be done in a forthcoming work. 
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APPENDIX A: PHOTON WAVE FUNCTION COMPUTATION 

The photon wave function in the framework of light cone perturbation theory is discussed 



m 



TB|. We compute it using the rules and conventions given in For a photon with 



momentum q = [q~^,q~ = —Q'^/2q~^, q = 0], one multiplies 

• a color factor ^/7^, 

• a flavor part ef6fj, 

• a spinor term u{zq'^, k, h) efj,{q, X)-y^v{{l — z)q'^, — k, h), 
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• a factor {\^zq^) ^/^(^/^(l — z)q~^) for the quark and antiquark lines, 

• a light cone energy denominator —\plq^ + ^^-^^ 

The polarization vectors of the photon are 0) = [g^/Q, Q/2g+, 0] and ±1) = 
[0, 0, — =fVv^]- The spinor matrix element between infinite-momentum-frame helicity 
eigenspinors p2[ are 



h,—hi 



/^(l - z)q+ 



- z) Jz{l - z) 



In the last line z = 1, 2 and =F stands for a minus sign if /i = +1/2 and for a plus sign when 
h = -1/2. 

By taking the Fourier transform 

one gets the expressions given in Eq. (p.2|) . In the longitudinal case there is an additional 
5'^^^(r) which one can drop because the color interaction vanishes at transverse distance. 

It is of course possible to obtain a wave function description in covariant approach and we 
want here to give the steps necessary to get the photon wave function. First, it is important 
to notice that a photon-quark-antiquark coupling in a Feynman graph can be interpreted 
in term of a photon wave function in light cone perturbation theory if the x"'"-ordering is 
7 ~^ QQ- This is the case in the formal limit — > +oo where this ordering survives. 

The first step is the evaluation of the k~ integral which leads for asymptotic q'^ to 

dk+d^k ^j^- f{k+,k-,k) dzd^k if{k+, 0-,k) 



(27r)^ j-oo {k"^ — + ie){{k — qY — rn? + is) IQn'^ k? + m? — z(l — z)q'^'' 

provided f{k~ = 0) is finite and non-zero. The numerator, = i{I^ + 'm){—ie^^)i{Iit~^ + m), 
is 

AT = ie {fc ■ e (2^ - ^) - ■ (fc - g) / + ie„^^p757°fc^e''g^ + m{2k ■ e - + m^)} 



6 



ieq+-f- I -^[z{l - z)q^ + + m"] + (1 - 2z)k ■ st + i^^e'^^k's^ + 



mi 



The wave function is obtained by taking the helicity matrix element w{h)Nw{—h) /\/2q~^ 
with w{l/2) = (1/V2, 0, -1/V2, 0) and ^(-1/2) = (0, 1/^2, 0, 1/^2) |2|. 
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APPENDIX B: HADRON WAVE FUNCTION PARAMETERS 



The value of the wave function at the origin is related to the meson leptonic decay 
constant 

(0|J'^(0)|F(g,A)) = e/vMv£'^(g,A), (Bl) 
which appears in the expression of the vector meson e~^e~ width 

In the parametrization Eq. ( |3.3| ) , this constraint leads to the determination of the parameter 
A/". The fixing condition, which depends on the meson helicity, A, is 

■rf..(l-.)W.)=/<;.^^^i§ffp^«.). (B2) 

Jo 2Myz{l — z) 

with ey the effective quark charge in the meson V expressed in units of the proton charge 
(see Table |). 

The normalization condition is 

{V{q\ nV{q. A)) = (27r)32g+5(g+ - e)5^'\^ - qO^aa', (B3) 
which leads to the relation 

^x = ^==:-yh, (B4) 



where 



Ji, = f^dzz^l-zffl(z), 



ujx and N\ are therefore defined by a system of implicite equations, Eq. ( P^] ) and Eq. 
Solutions are lis ted in Table H together with the corresponding root mean square radius 
R= (i?2)i/2 = J^l2uj. 
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